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■ ABSTRACT 

o : 

^ . We present mid-infrared ISO-SWS observations of the fine-structure emissions 

'. lines [Ne II]12.8/im, [Ne III]15.6;um, [Ne III]36.0/im, [Ar II]6.99/im, [Ar III]8.99/im, 

> . [S III]18.7/im, [S III]33.5yum, [S IV]10.5/im and the recombination lines Bra 

^ . and Br/3, in a sample of 112 Galactic H II regions and 37 nearby extra- 

galactic H II regions in the LMC, SMC, and M33. We selected our sources 
from archival ISO-SWS data as those showing prominent [Ne II]12.8/im or 
[Ne III]15.6/im emissions. The Galactic sources have a wide range in galacto- 
centric distance (0 kpc <, Rgai ^ 18 kpc), which enables us to study excita- 
tion and metallicity variations over large Galactic scales. We detect a steep 
rise in the [Ne III]15.6/im/[Ne Il]12.8/im, [Ar IIl]8.99/im/[Ar Il]6.99/xm, and 
[S IV]10.5/im/[S IIl]33.5/im excitation ratios from the inner Galaxy outwards, 
and a moderate decrease in metallicity, from ~ 2Zq in the inner Galaxy to 
~ 1^0 in the outer disk. The extra-galactic sources in our sample show low gas 
density, low metallicity and high excitation. We find a good correlation between 
[Ne III]15.6/im/[Ne II]12.8yum and [Ar IIl]8.99/im/[Ar II]6.99/im excitation ratios 
in our sample. The observed correlation is well-reproduced by theoretical neb- 
ular models which incorporate new generation wind-driven NLTE model stellar 
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atmospheres for the photoionizing stars. In particular, the NLTE atmospheres 
can account for the production of [Ne III] emission in the H II regions. We have 
computed self-consistent nebular and stellar atmosphere models for a range of 
metallicities (0.5 to 2 Zq). We conclude that the increase in nebular excitation 
with galactocentric radius is due to an increase in stellar effective temperature 
(as opposed to a hardening of the stellar SEDs due to the metallicity gradient). 
We estimate an integrated [Ne III]15.6/im/[Ne II]12.8/xm ratio for the Galaxy 
of 0.8, which puts it well inside the range of values for starburst galaxies. The 
good fit between observations and our models support the conclusion of Thorn- 
ley et al. (2000), that the low [Ne IIl]15.6/xm/[Ne Il]12.8//m ratios observed in 
extra-galactic sources are due to global aging effects. 

Subject headings: galaxy: abundances — galaxy: stellar content — nebulae: H II 
regions — stars: atmospheres 



1. Introduction 

Mid-infrared (IR) fine-structure emission lines are useful probes of the physical proper- 
ties of H II regions and their associated ionizing sources. Unhke optical and UV emission 
lines, the fine-structure lines do not suffer as much from dust extinction, and the low en- 
ergies of the fine-structure levels {hu/k ~ 10^~^K), make the fine-structure hnes relatively 
insensitive to the nebular electron temperature Te (~ lO^K in H II regions). Fine-structure 
line emissions from photoionized gas have been observed in a wide range of astrophysical 
sources (Genzel et al. 1998; Cox ct al. 1999; Thornley et al. 2000). Their measured 
intensities have been shown to be strong, comparable to optical hydrogen recombination 
lines. The relative strengths of the fine-structure emission lines may be used to probe the 
ionizing spectral energy distributions (SEDs) because different ionization stages are formed 
by significantly different photon energy ranges in the ionizing spectrum. For example, the 
[Ne IIl]15.6/im/[Ne Il]12.8//m line ratio is sensitive to photons emitted at > 3 Rydberg re- 
gion, while [Ar IIl]8.99//m/[Ar Il]6.99/im probes the region > 2 Rydberg. Measuring these 
(and other) ratios can constrain the SEDs of the hot stars in H II regions (Thornley et al. 
2000; Forster Schreibcr et al. 2001a) or of the non-stellar sources in active galactic nuclei 
(Alexander et al. 1999; 2000). 



^Based on observations with ISO, an ESA project with instruments funded by ESA Member States (espe- 
cially the PI countries: Prance, Germany, the Netherlands, and the United Kingdom) with the participation 
of ISAS and NASA. 
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The Short Wavelength Spectrometer (SWS) on board Infrared Space Observatory (ISO) 
made spectroscopic studies of Galactic and extragalactic sources possible across the entire 
mid-infrared (2.5-45/im) spectral range (de Graauw et al. 1996; Kesslcr ct al. 1996; Genzel & 
Cesarsky 2000). In this paper we have compiled SWS observations of fine-structure emission 
lines from Galactic and nearby extra-galactic H II regions from the ISO data archive. The 
sources were selected as those in which the [Ne II]12.8/xm or [Ne III]15.6/xm lines are present, 
and eliminating target types other than H II regions. Since these objects are typically 
prominent sources of the diagnostically useful [Ar II]6.99//m, [Ar III]8.99//m, [S III]18.7//m, 
[S IIl]33.5/xm, and [S IV]10.5//m fine-structure fines and the Bra and Br/3 recombination 
lines, we extracted these lines as well for sources in which they were observed. 

In this paper we have several goals. First, we wish to re-investigate the large scale 
variations in elemental abundances and nebular excitation across the Galactic disk. The 
Galactic sources in our sample have a wide enough range in galactocentric distance (0 kpc ^ 
Rgai ^ 18 kpc) to enable us to study these variations. Such studies using ISO data have 
been presented by Cox et al. (1999) for a sample of 46 Galactic H II regions over a 
range of kpc to 25 kpc in galactocentric radius, and showed that the excitation ratios 
[Ne IIl]15.6/Lim/[Ne Il]12.8/im and [Ar IIl]8.99;Lim/[Ar Il]6.99/im increase as a function of 
the galactocentric radius. They also showed that the argon and neon abundances, based on 
observations of [Ne II]12.8/xm and [Ar II]6.99/i.m decrease by a factor of 2 from the inner 
Galaxy outwards (see also Roelfsema et al. 1999 and Martm- Hernandez 2000). These results 
are consistent with the trends found in earlier studies of Galactic H II regions using ground- 
based and air-borne observations (Simpson, Colgan, & Rubin 1995; Afflerbach, Churchwell, 
& Werner 1997). Here we present results for a larger source sample. Understanding the ex- 
citation gradient within our galaxy may shed further light on excitation variations observed 
in external galaxies (e.g. Shields & Tinsley 1976 in MlOl; Henry & Howard 1995 in M33, 
M81 and MlOl; see also Shields 1990). 

Second, we use the observations as templates for constraining the results of new gen- 
eration wind- driven (NLTE) stellar atmosphere computations (Pauldrach, Hoffmann, & 
Lennon 2001). For example, a long-standing problem in Galactic H II regions studies has 
been the failure of photoionization simulations employing old-generation static (LTE) at- 
mospheres (Kurucz 1992) to reproduce the observed (high-excitation) [Ne III] line emission 
(Baldwin et al. 1991; Simpson et al. 1995). We re-investigate (and resolve) this issue here 
using NLTE atmospheres. In our analysis we construct an [Ar III]8.99//m/[Ar Il]6.99//m vs. 
[Ne III]15.6//m/[Ne II]12.8/im diagnostic diagram as a measure of the nebular excitation. We 
fit the observed correlation between these two line ratios by constructing nebular simulations 
in which the photoionizing stellar SEDs are modeled using NLTE atmosphere computations. 
We compare with results obtained using the older LTE Kurucz atmospheres. We also explore 
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the effects of metallicity variations by self-consistently varying the nebular and the stellar 
atmosphere abundances. 

This paper is a part of a larger framework, whose aims are to study theoretical models 
of starburst regions, as they are probed by IR spectroscopy (Kunze et al. 1996; Lutz et 
al. 1996; Rigopoulou et al. 1996; Sturm et al. 1997; Genzel et al. 1998; Thornley et al. 

2000, Forester Schreiber et al. 2001a. b). For example. Thornley et al. (2000) measured the 
[Ne IIl]15.6/im/[Ne Il]12.8/im excitation ratio in extra-galactic starburst sources, and found it 
to be low compared to starburst models which incorporate the new-generation atmospheres. 
They concluded that this is evidence for global aging of the star-forming regions. Fine- 
structure emission lines have also been used to try and determine the primary energy source 
in the ultra-luminous infrared galaxies (ULIRGs; Lutz et al. 1996; Genzel et al. 1998; 
Rigopoulou et al. 1999). Understanding the properties of the SEDs of hot stars, as probed 
by Galactic H II regions, is crucial to the interpretation of the extra-galactic observations. 



2. Observations 

2.1. Data Reduction and Source Selection 

The ISO data archive contains several thousands of observations that were obtained 
with the 2.4-45/im grating spectrometer SWS. We have identified potential H II regions in 
all observations taken in the full spectrum mode (SWSOl) and the line spectroscopy mode 
(SWS02) during the ISO mission. An automatic procedure starting from the off-line pro- 
cessing version 7 standard processed data (SPD) processed all these observations to the auto 
analysis result (AAR) level. From every AAR, a spectral range of ±30 resolution elements 
(r=2000) around the [Ne II]12.8/im and [Ne III]15.6/im hues, was cut out for further analysis. 
This processing was automated by a script which flat-fielded the individual detectors to their 
mean value, removed 3a outliers and rebinned each line to a single- valued spectrum. 

For the SWSOl data, a simple automatic removal of the SWS band-3 instrumental 
fringes, was applied in this procedure when the continuum for the [Ne II]12.8/xm and 
[Ne IIl]15.6//m fines exceeded 100 Jy. We then automatically fitted gaussians to the re- 
binned spectral segments and measured line fiuxes. A source was considered a candidate for 
in-depth processing if at least one of the two lines was apparently present at a level above 
10-20 Wcm-2. 

For the SWS02 data, the automatically produced single valued spectral segments for 
each line were stored and later inspected visually for the presence of lines, which finally gave 
a list of observation number and target name for every clear detection of either [Ne Il]12.8//m 
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or [Ne IIl]15.6//m or both. Visual inspection was preferred here, since the defringing step in 
AOT02 observations is less robust, and a lot of faulty detections would occur otherwise. 

Our selection procedure initially gave also objects of similar spectral characteristics 
(nearby galaxies, planetary nebulae etc.) which were subsequently identified and excluded 
from further processing. All observations from the hst of candidate H II regions were then re- 
processed and the line fluxes were measured. In Table 1 we list the the specific emission lines 
included in our source sample, together with their ionization potentials, line wavelengths, 
and critical densities for electron impact collisions. 

For the SWSOl data, this needed special care since a large fraction of the observa- 
tions was taken in the faster speed options of the SWSOl mode. For those observations, 
the grating moved during integrations causing curved integration ramps when passing a 
strong emission line. The automatic deglitching of the SWS data occasionally confused 
such curved ramps with jumps in the ramps that correspond to glitches caused by ion- 
izing particle hits, and thus might have affected emission hne fluxes (Section 4.5 of ISO 
Handbook Volume 6, The Short Wavelength Spectrometer, SAI-2000-008/Dc, Version 1.1., 
http://www.iso.vilspa.esa.es/manuals/HANDBOOK/VI/sws_hb/). All SWSOl data except 
for the slowest speed 4 were hence reprocessed twice from ERD level, once with standard 
OLP Version 8 calibration files and once with one calibration file changed to switch off the 
glitch recognition. In both cases, the data were processed to the AAR level, cleaned from 
outliers by iterative clipping, flat-fielded, and rebinncd to a single-valued spectrum. The 
two versions were then plotted and inspected to determine whether the deglitching might 
have affected the real line in the standard processing, or whether a glitch might have been 
present in those data with glitch recognition switched off. This classification was in most 
cases possible on the basis of feature shape and widths. The line fiux or limit was then mea- 
sured in that version not affected by problems. In a few cases with both versions affected, 
generous upper limits were adopted. SWSOl speed 4 and SWS02 data were processed by 
default before measuring line fiuxes by direct integration of the continuum subtracted line, 
or determining upper limits corresponding to a line, of 3(J peak height and the proper width. 
For SWSOl observations this was done for all lines in Table 1, for SWS02 only for those 
present in a particular observation. The uncertainty in the line fiuxes is estimated to be 
~ 30%. 

The final identification of individual objects - whether it is an H II region or a different 
kind of source (e.g. nearby galaxy, planetary nebula etc.) that was missed by the flrst 
filtering, or whether there were multiple observations of the same object - is based on a 



-6- 



literature survey, mainly using the SIMBAD astronomical database^. 

A total of 149 SWS spectra were selected for the final sample. This includes almost all 
of the H II regions observed by ISO-SWS - 112 Galactic H II regions, spanning galactocentric 
radii of to ~18 kpc, 14 H II regions in the Large Magellanic Cloud (LMC), 4 in the Small 
Magellanic Cloud (SMC), and 19 in the nearby spiral galaxy M33. This sample is certainly 
heterogeneous since it is a collection of targets observed for different proposals and scientific 
reasons. These include detailed studies of well known individual H II regions, more systematic 
surveys of compact H II regions originally selected on the basis of IRAS-LRS spectra and 
on the basis of IRAS sources located outside the solar circle, observations of the brightest 
H II regions in the LMC, SMC, and M33, and H II region observations targeted at specific 
issues of ISM physics rather than a per se characterization of the H II region observed. The 
selection criteria adopted by these proposals necessarily differed, and incomplete execution of 
planned observations during the ISO mission induces another (though uncritical) selection 
effect. Some observations represent integrated fluxes of rather compact well-defined H II 
regions, others different parts of more complex regions, including pointings in the outskirts 
of resolved nebulae. However, the sample spans H II regions with a wide enough range of 
positions and excitation conditions to study both the stellar SEDs and Galactic metallicity 
and nebular excitation gradients. 

Metallicity gradient analyses in our Galaxy should be little affected by possible selection 
biases in our sample, as long as there is efficient mixing at a given galactocentric radius. The 
situation is less clear for excitation analyses where competing selection biases are possible. 
Selection on compactness, for example, might prefer younger (harder) but also less luminous 
(softer) objects, with additional effects due to dust possible. Similarly, optical rather than 
infrared/radio selection might prefer more evolved (softer) but also more luminous (harder) 
and lower metallicity (harder) objects, all perhaps more likely to be optically visible. Since a 
large fraction of the targets in the ISO database is originally IRAS-selected, such effects are 
likely not governing excitation analyses as a function of galactocentric radius in our sample. 

2.1.1. Calculation of Galactocentric Radius 

We obtained radial velocity measurements for most of the Galactic sources from the 
hterature (see references in Tables 2 and 3). These velocities were determined using dif- 
ferent methods and are based on various emission features (e.g., radio recombination lines; 
sub-mm molecular hues of CO, CS, and CI; maser emission). In some cases (DR21, GAL 



^ http : / /simbad . u-str asbg . fr / 
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007.47+00.06), the given velocity value is the mean between different values obtained from 
the same method (i.e., several recombination lines). 

We use the Rohlfs & Kreitschmann (1987) Galactic rotation curve to calculate the 
galactocentric distance of each source. Their rotation curve spans a range of 0.1 to 19.4 
kpc in galactocentric radius. For the Sun we adopt a galactocentric distance Rq — 8.5 kpc 
and rotational velocity 6q = 220 km s~^ (Honma & Sofue 1997). We have transformed the 
rotation curve of Rohlfs & Kreitschmann {Rq = 7.9 kpc and 6q = 184 km s^^) using the 
more updated solar values of Honma & Sofue. Most of the Galactic sources lie very near the 
Galactic plane (typically, 6 < 5°), so we assume that they lie in the Galactic plane (6 = 0). 
(An exception is the Orion nebula (6 = —19.4"), for which the following procedure was not 
applied.) For each source having a radial velocity vios and a Galactic longitude I, we compute 
the quantity W{Rgai), where vios — W{Rgai) sinZ, and 



Here, Rgai is the galactocentric radius and 9{Rgai) is the rotational velocity at Rgai- Using 
a linear interpolation of the W{Rgai) curve from Honma & Sofue we find the values of Rgai 
corresponding to the W{Rgai) of each source. 

The kinematic distance {Rkm, distance from the Sun) is calculated by inverting the 
formula 



For sources outside the solar circle, this gives a single solution, whereas inside the solar circle, 
one obtains two solutions. The choice of distance in this case is based on different methods 
discussed in the references. A distance ambiguity remains for only three sources designated 
with asterisks in Table 2 (column 8). 



Table 2 lists the data collected from the hterature on the Galactic sources. The columns 
in this table give the following information: 

Column ( 1 ) — Source identification number in the sample. These numbers also appear 
as source identifications in the various figures. 

Column (2) — Source name. 




(1) 



R^al — "^RkinRe COS / + 



(2) 



2.1.2. The Sources Data 
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Column (3-4) — Position (J2000) of the ISO aperture in Right Ascension and dechna- 
tion. 

Column (5) — Kinematic distance (distance from the Sun) in kpc. 
Column (6) — Galactocentric radius (see § 2.1.1) in kpc. 
Column (7) — Source radial velocity in km s~^. 

Column (8) — Acronym for the distance/radial velocity reference. The literature refer- 
ences are given in Table 3 

Column (9-18) — ISO-SWS fluxes for the lines listed in Table 1. The fluxes are in units 
of 10~^^ W cm~^. The error in the flux measurement is ~ 30%. 

In several cases, where observations of different parts of spatially resolved and complex 
regions were available, we summed the multiply observed intensities across the regions to a 
single value for each line. These observations are listed in Table 4. Table 5 hsts the data for 
sources in the LMC and SMC. In Table 6 we list the data for sources in M33. In this table 
we also list the angular distance of each source from the center of the M33. 

A few studies of individual sources in our sample have been presented in the litera- 
ture. For these sources our fluxes are in good agreement with the previously published 
fluxes. For example, in Orion IRC2 van Dishoeck et al. (1998) give fluxes for all of our 
lines (except for Br/3) which are similar to the fluxes of our source #13. Similarly, van den 
Ancker, Tielens, & Wesselius (2000) pubhshed fluxes for the [Ne Il]12.8/xm, [Ne IIl]15.6//m, 
[S IIl]18.7/im, [S IIl]33.5;um, and [S IV]10.5yum fine-structure lines of S106 IR, which are com- 
parable to those of our S106 1RS4 (#99). Rosenthal, Bertoldi, & Drapatz (2000) published 
fine-structure fluxes for all of the lines we use in this paper for the Orion H2 emission peak 
nebula (our sources #11 and #12). Rodriguez-Fernandez, Mart in- Pint ado, & de Vicente 
(2001) give fluxes of several sources in the Galactic center from which six sources are included 
in our sample (#62-#67). They measured [Ne II]12.8/im, [Ne III]15.6/im, [S III]18.7/im, and 
[S III]33.5//m. 

3. Results 
3.1. Excitation 

We use the line ratios [Ne IIl]15.6/xm/[Ne Il]12.8//m, [Ar IIl]8.99/xm/[Ar Il]6.99//m, 
and [S IV]10.5//m/[S III]33.5//m to quantify the source excitation. In Figure 1 we show 
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histograms of the observed [Ne IIl]15.6//m/[Ne Il]12.8//m ratio (not corrected for extinction) 
for the groups of sources in our sample - the Milky Way, the LMC, the SMC, and M33, 
and for comparison, the sample of extra-galactic sources, as measured by Thornley et al. 
(2000). Filled bars designate sources whose ratios are uncertain. For the Milky- Way and 
the starburst sample these are all upper limits, for the LMC and SMC they are lower limits, 
whereas in M33 both cases are mixed together. It appears that the excitation distribution 
of Galactic H II regions is quite different from the distribution in the LMC, SMC and M33, 
with smaller excitation in the Milky Way on average. 

We can use our observations of Galactic H II regions to estimate an integrated value 
of the [Ne III]15.6/xm/[Ne II]12.8/im ratio for our Galaxy which can be compared to the 
starburst galaxies. This is necessarily approximate since the sample drawn from the ISO 
archive is not an unbiased representation of the nebular luminosity of the Galaxy. The 
integrated ratio (not corrected for extinction) equals 0.8, with a de-reddened value (see §3.2) 
of 0.7. The corresponding median values are 0.41 and 0.36. These, admittedly approximate 
values (shown as vertical lines over-plotted on the Galactic and extra-galactic histograms in 
Figure 1) put our Galaxy well inside the range of [Ne III]15.6//m/[Ne II]12.8/im values for 
much more active starburst galaxies observed by Thornley et al. (2000), and even above its 
median. 

In Figure 2 we plot the three excitation ratios mentioned above vs. galactocentric 
radius for the Galactic sources. The data points are represented by the source identification 
numbers, as listed in Table 2. The displayed ratios have been corrected for extinction (§ 3.2). 
In this figure and the figures that follow we designate limits on the given quantities by small 
arrows. Large arrows designate lower limits on the extinction correction. 

Figure 2 shows a clear rise in the excitation towards the outer Galaxy. We have fitted 
power- law curves to the data in the range 2 < Rg^i < 11 kpc. The fitted gradients are 

and 

The observed excitation gradients are similar to those found in earlier studies of Galactic 
H II regions (Simpson et al. 1995; Afflerbach et al. 1997, Cox et al. 1999). We note that the 
[Ne IIl]15.6/xm/[Ne Il]12.8//m and [Ar IIl]8.99/xm/[Ar Il]6.99//m gradients are quite similar. 
This is an important clue as to the origin of the excitation gradients as we discuss in § 4. 
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The ionization potentials of Ar+, S++, and Ne+ are equal to 27.63 eV (~ 2 Ryd- 
berg), 34.79, and 40.96 cV (~ 3 Rydbcrg), respectively (see Figure 6). Thus, the ratios 
[Ar III]8.99/im/[ArII]6.99/xm, [S IV]10.5/im/[S III]33.5/im, and [Ne III]15.6/xm/[Ne II]12.8/xm 
constrain the shapes of the SEDs of the ionizing spectra over a range of several Rydbergs. 
One possible representation for the data, which may give clues regarding the properties of the 
SEDs, is a diagnostic diagram combining two of these empirical ratios. In Figure 3 we plot 
the [Ar III]8.99/im/[Ar II]6.99/im vs. [Ne III]15.6/im/[Ne II]12.8/xm ratios for our sample. A 
good correlation is apparent, showing that nebular excitation is probed by both line ratios, 
and manifests the inter-relation between the two spectral regions, 2 < £ < 3 Rydberg and 
£ > 3 Rydberg. A representative error bar for the data is shown in the lower right corner of 
the plot. 

In order to connect this correlation in excitation space to position within the Galaxy we 
display each source in the lower panel of Figure 3 as a symbol according to its galactocentric 
radius. Pluses designate the inner kpc in the Galaxy {Rgai < 1 kpc), circles designate 
Galactic disk sources {Rgai > 1 kpc), and asterisks designate extra-galactic sources. From 
this diagram it is clear that the sources in the central kpc of the Galaxy, are all of low 
excitation. On the other hand, sources in the LMC, are all of high excitation. Sources in 
the disk of the Galaxy show moderate levels of excitation. 

3.2. Extinction and Gas Densities 

We use the Brackett line ratio, Bra/Br/S (for those sources for which both lines are 
measured) to de-redden the fine- structure line intensities. We assume a foreground screen, 
and adopt the Draine (1989) mid-IR extinction curve. For case-B recombination the intrinsic 
Brackett decrement ranges from 1.6 to 1.9 for gas temperatures between 5 x 10^ and 2 x 10^ 
K, and densities between 10^ and 10^ cm~^. Figure 4 shows the Brackett ratio for the objects 
in our sample vs. galactocentric radius. 

It is evident that there are objects with unphysical Brackett ratios much below the 
optically thin ratio of ~ 1.7. These objects generally correspond to sources with noisy low 
signal-to-noise spectra, and we discard them in our subsequent analysis. The extinction 
correction is most significant for the [Ar III]8.99//m and [S IV]10.5//m fines which lie close 
to the mid-IR silicate absorption feature. For these fines correction factors as large as 10 
to 25 are implied by the Brackett decrements. The corrections for the other fine-structure 
lines are small. For some objects only an upper-limit exists for Br/3. We adopt the upper- 
limit values to compute a minimum extinction correction for such objects. Objects with 
lower-limit extinction corrections are indicated by the large arrows in the figures. 
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The [S IIl]18.7//m/[S IIl]33.5/xm ratio was used as a gas density measure. The sohd 
angle of the SWS aperture used to observe the 33.5/xm hne is 1.75 times larger than the 
one used for the 18.7/im line. Depending on source structure, this may require aperture 
corrections to the ratio up to about this factor. We have not attempted to derive these 
corrections for individual sources but note that observed ratios of ftiO.25 are possible for low 
density extended sources. 

In Figure 4 wc display the [S IIl]18.7/im/[S IIl]33.5/um vs. galactoccntric radius for the 
(extinction corrected) Galactic sources. The measured densities for our sample lie in the 
range of <^ 100 to 15,000 cm~^ with a mean value of 800 cm~^. Measurements of the 
S III ratio are also available for some of the LMC sources. These ratios lie in the range of 
0.4 — 0.6, corresponding to densities ^ 200 cm"^. 



3.3. Abundances 

We have calculated the elemental abundances for neon and argon from the relative 
strengths of the fine-structure and Bra emission lines. In computing the fine-structure line 
emissivities we adopt the "IRON project" collision strengths (Saraph & TuUy 1994 for [Ne II]; 
Butler & Zeippen 1994 for [Ne 111]; Pelan & Berrington 1995 for [Ar II]; Galavis et al. 1995 for 
[Ar III]). (These data are also employed in the CLOUDY photoionization models we discuss 
in §4.) The critical densities (> 10^ cm~^) of the relevant neon and argon transitions are 
much larger than the estimated gas densities so that coUisional de-excitations are negligible. 
Because of the large ionization potential (4.7 Ryd) of Ne"*"^ we assume that neon exists 
either as Ne"'" or Ne"'""'" in the H II regions. The ionization potential of Ar"*"^ is 3.0 Ryd, 
very similar to that of Ne+. However, Ar "'"+"'" ions are rapidly removed by charge transfer 
reactions with neutral hydrogen (Butler & Dalgarno 1980). In contrast, charge transfer 
neutrahzation is inefficient for Ne"'""'" (Butler, Heil, &: Dalgarno 1980). We therefore assume 
that negligible [Ar IV] is produced in the nebulae even when [Ne lll]/[Ne 11] is large. This is 
verified by our CLOUDY computations which incorporate the Butler et al. charge transfer 
rate coefficients. Thus, we do not apply ionization "correction factors" in deriving the neon 
and argon abundances from our fine-structure line data. 

For the singly ionized species, Ar"'" and Ne"'", which have two fine-structure levels in their 
ground state, the emissivity is given by 

jii = nenqiihuii^, (6) 
rl 

where n is the gas density,ne is the electron density, qu is the coUisional excitation rate, and 
Xn/H is the abundance of the species relative to hydrogen. For the doubly ionized species. 
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Ar+''" and Ne++, which have three fine-structure levels in their ground state, the emissivity 
is given by 



Jill = n^n{qi2 + ^isjAii^iii-rr = n^nqinhvin- 



(7) 



H "'^"■^^^ H ' 

where qyi and ^13 are the coUisional excitation rates between the first level and the second 
and third levels, respectively. The Bra emissivity 



(8) 



where nj^+ is ionized hydrogen density 0:^(5 4) is the case B effective recombination 
coefficient for Bra. We assume n ^ nH+- 



It follows that the abundance, X/H, is given by 
X ^ Xn _^ Xni ^ a(5 4)/it/Bm 

H H H jBra 



+ 



JlII 



qiihuii qiiihuiii^ 



(9) 



Inserting the basic atomic data assuming a gas temperature of lO'^K one gets: 



Ne 

IT 



j^Q_5 ^ 2.30/Ne II + /Ne III 



1.90/, 



Bra 



Ar 

IT 



j^Q_6 ^ /Ar II + 1-26/Ar III 



1.43/, 



(10) 



Bra 



where fx are the line ffuxes. 



Figure 5 shows the computed neon and argon abundances vs. galactocentric radius for 
the Galactic sources. Both abundances show a decrease from the inner parts of the Galaxy 
outwards. Both elements vary from ~ 2x solar to ~ Ix solar metallicities. The fitted 
gradients over the entire range of galactocentric radii are 



and 



Ne' 

IT 



Ar 

IT 



where the solar metaUicity values are [Nc/H]q 
Sauval 1998). 



= -3.71(±0.06) - 0.021(±0.007)i?g„,, 

= -5.31(±0.06) -0.018(±0.008)i?g„,, 
-3.92 and [Ar/H]o 



(11) 

(12) 

-5.60 (Grevesse & 



The statistical significance of our detected gradients is low. However, our results are 
consistent with previous studies. Afflerbach et al. (1997) found similar gradients, from twice 
solar metaUicity in the inner Galaxy to lower values outwards, in their study of 0/H and 
S/H. Simpson et al. (1995) found a gradient in the neon abundance which is steeper than 
ours (— 0.08±0.02 dex kpc~^), but their sample consists of only 18 sources. Simpson & Rubin 
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(1990) also found the inner Galactic sources to be of about twice solar metallicity. In their 
preliminary study Cox et al. (1999) show evidence of a neon abundance gradient as probed 
by the [Ne II]12.8/im/BrQ; ratio. However, their analysis did not include contributions due to 
Ne"*""*" ions as probed by the [Ne III] emission which dominates in the high excitation nebulae. 

4. Analysis 

The aim of the following analysis is to obtain a theoretical fit to the empirical 
[Ar IIl]8.99/xm/[Ar Il]6.99/xm vs. [Nc IIl]15.6/im/[Ne Il]12.8/xm diagram shown in Figure 3. 
We consider a sequence of photo-ionization models in which we assume that the H II re- 
gions are ionized by single stars. We compute NLTE model atmospheres for high surface 
gravity dwarfs and low gravity giants stars with a range of effective temperatures. In the 
nebular analysis we explore the effects of varying the gas density, ionization parameter and 
metaUicity. We also consider the effects of metaUicity variations in the model atmosphere 
computations. 

4.1. Models 

We use the WM-basic stellar atmosphere codc"^ (Pauldrach ct al. 2001) to compute the 
spectral energy distributions. These NLTE atmosphere computations include the hydrody- 
namical treatment of steady-state winds and mass loss, and include the combined effects of 
line blocking and line blanketing. The WM-basic package employs up-to-date atomic data, 
which provide the basis for a detailed multi-level NLTE treatment of the metal ions, from C 
to Zn. As discussed by Pauldrach et al., a key (and unique) feature of our model atmospheres 
is that they successfully reproduce the spectral characteristics of hot stars as revealed by 
high-resolution far-UV spectroscopy. We use the NLTE atmospheres to construct model 
stars as input ionizing sources for the nebulae which we model using the photo-ionization 
code CLOUDY (version C90.04, Ferland et al. 1998)*^. We also carry out computations using 
the older generation LTE static atmospheres (Kurucz 1979, Kurucz 1992). 

In general, the NLTE wind models have "harder" SEDs compared to the hydrostatic 
LTE models (see also Schaerer & de Koter 1997). For a given stellar effective temperature, 
the NLTE atmospheres emit significantly larger fluxes of energetic photons. Figure 6 shows 



"^Sec also www. usm.uni-muenchcn.dc/pcople/adi/adi. html 
''See also www.pa.uky.edu/~gary/cloudy/ 
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an example of the differences between the LTE and NLTE spectral energy distributions for 
an effective temperature Tg// = 35, OOOK and surface gravity log 51 = 4.0. 

The relative deficiency in high energy photons 3 Rydberg) in the LTE model 
is clear. The ionization edges of the relevant ions are overlayed on the spectra in this 
plot. The nebular ionization state, expressed by the ratios [Ne III]15.6//m/[Ne II]12.8//m or 
[Ar IIl]8.99/xm/[Ar Il]6.99//m is sensitive to the shape of the SED, and the NLTE models 
will generally produce greater excitation. It is apparent from Figure 6, that the difference 
in excitation produced by the LTE and NLTE atmospheres will be more pronounced in the 
Ne"'"+/Nc+ ratio compared to Ar++/Ar+. This helps in breaking the degeneracy between 
the SED shape and the ionization parameter as factors which influence the ionization state 
of the gas. 

We have computed model atmospheres for "representative" main sequence dwarfs 
(loggf = 4.0) as well as super-giants {logg = 3.25 — 3.75) for T^fj ranging from 35,000 
to 45,000 K, in steps of 10^ K. We computed models for solar (1 Zq; Grevesse & Sauval 
1998), half solar (0.5 Zq), and twice solar (2 Zq) metallicities. The stellar wind opacities 
and resulting SEDs depend on the assumed stellar mass loss rates (Pauldrach et al. 2001). 
For each metallicity we fix the mass-loss rates (and terminal wind velocities) using the em- 
pirically based wind momentum-luminosity relationship for 0-star super-giants, giants, and 
dwarfs (Kudritzki & Puis 2000). 

We represent the nebulae as optically thick Stromgren spheres, and consider clouds 
with constant hydrogen densities ranging from 10^ to 3 x 10^ cm~^. We assume a gas filling 
factor f — 1. For each stellar effective temperature, Tg//, we compute the Lyman continuum 
photon emission rate, 

Q ^ A^Rl ^^^^^^dv is-'), (13) 

where Fy{T(.ff) (in erg s~^ cm ~^ Hz~^) is the SED, i/q is the Lyman limit frequency, and i?* 
is the stellar radius (~ 10 Rq for 0-typc dwarfs, and ~ 20 Rq for giants). 

The nebular ionization structure depends on the shape of the SED, and also on the 
ionization parameter 

where c is the speed of light, n is the nebular gas density, and is the Stromgren radius 
given by 

47r 

Q = yasr^V, (15) 
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where as is the case B recombination coefficient. It follows that 

Thus, for a fixed nebular density the ionization parameter is fully determined by the stellar 
atmosphere model. For a given density we determine U from the stellar model, and compute 
the resulting [Ar III]8.99/im/[Ar II]6.99/im and [Ne III]15.6/im/[Ne II]12.8/xm emission line 
ratios using CLOUDY. As a representative density we set n — 800 cm~^ as determined from 
the S III density diagnostic (see §3.3). 

Tables 7, 8, and 9 summarize the parameters for the stellar atmospheres we adopt as the 
ionizing sources in our photoionization computations. These parameters include the stellar 
effective temperature, surface gravity, stellar radius, wind parameters (mass-loss rate M and 
the terminal velocity v^), and hydrogen photoionizion rates for NLTE typical dwarfs and 
super-giants, and for LTE dwarfs. The parameters are given for the half-solar, solar, and 
twice-solar metallicities. 



4.2. Comparison with Observations 

In Figure 7 we show the [Ar IIl]8.99//m/[Ar Il]6.99//m vs. [Ne III]15.6//m/[Ne Il]12.8/im 
observations (Figure 3) together with the model computations for the LTE (dwarfs) and 
NLTE (dwarfs and giants) stars for solar metallicity, and a gas density of 800 cm"^. The 
bold faced numbers beside each model point give the T^ff of the corresponding star in units 
of 10^ K. A representative error bar for the data is shown in the lower right corner of the 
plot. 

For all sequences the predicted nebular excitation increases with increasing effective 
temperature of the photoionizing star. However, the LTE and NLTE models are clearly offset 
from one another, and it is evident that the NLTE sequences provide a much better fit to 
the observed excitation correlation. For a given [Ar III]8.99/xm/[Ar Il]6.99/xm ratio the LTE 
models predict [Ne III]15.6/im/[Ne II]12.8/xm ratios which are factors of ~ 10 smaller than 
observed. This discrepancy appears to be almost entirely ehminated in the NLTE sequences. 
The NLTE fit may actually be even better than is indicated by Figure 7 since some of the 
scatter is hkely due to underestimated extinction corrections for the [Ar IIl]8.99//m fine. 

The key difference between the LTE and NLTE models lies in the different spectral 
shapes of their photoionizing continua. This is illustrated in Figure 8 where for each ef- 
fective temperature in the LTE and NLTE sequences we plot the fractions Qargon/Q vs. 
Qneon/Q, whcrc Q argon IS the photou cmission rate for photons above the Ar+ 27.6 eV ion- 
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ization threshold, Qneon is the emission rate above the Ne+ 40.96 eV threshold, and Q is 
the total Lyman continuum (Lyc) emission rate. The SEDs harden as T^ff increases, and 
both ratios Q argon/ Q and Qneon/ Q iucrease. For a given value of Q argon/ Q the LTE at- 
mosphere models produce much smaller values of Qneon/ Q than the NLTE atmospheres. 
This is the behavior that is essentially represented in the [Ar IIl]8.99/im/[Ar II]6.99/Lim vs. 
[Ne III]15.6/xm/[Ne II]12.8/xm excitation diagram. The ISO fine-structure data thus provide 
direct empirical evidence for the improved accuracy of the NLTE atmosphere SEDs in the 
Lyman continuum range. The NLTE sequence also indicates that for most of the displayed 
H II regions the effective temperatures of the exciting stars plausibly lie in the range of 
35,000 to 45,000 K. 

We have examined the dependence of the model results on the ionization parameter by 
considering variations in both Q and n (see eq. 16). Although Q is fixed by the stellar model, 
we have verified that changing Q (for a fixed SED) simply moves the computed points along 
the (approximate) diagonal lines defined by each sequence of models. Larger values of Q 
could be produced by clusters of (identical) stars of given effective temperature. For an H II 
region excited by a mix of stars one would expect the high SED energies to be dominated 
by hotter stars and the low energies by cooler stars. We will explore the efi^ects of realistic 
clusters containing a synthesized distribution of stars and associated cluster SEDs in a later 
paper. Variations in the assumed density (in the range of 10^ to 3 x 10^ cm~^) induce only 
small shifts in the positions of the computed model sequences. 

In Figure 9 wc show the NLTE dwarfs and LTE model sequences for nebular densities 
of 100, 800, and 3 x 10^ cm^^. We conclude that the LTE sequence cannot be reconciled 
with the data by varying either Q or n. 

In order to investigate the influence of varying the metallicity we show in Figure 10 
results for three sequences of NLTE atmosphere models of dwarfs - computed assuming 
half-solar, solar, and twice-solar metallicities. We varied the metallicity self-consistently in 
the nebulae and in the stellar atmospheres. Both of these variations affects the computed fine 
ratios. Lower nebular metallicity results in higher gas temperatures, reduced recombination 
rates, and therefore slightly higher excitation. Afflerbach et al. (1997) found a similar 
trend in their study, but they also stressed the need for exploring the effects of metallicity 
variations in the stellar atmospheres as well. We have done this here, and find that such 
variations have a somewhat larger effect on the computed line ratios. For a given effective 
temperature the SED hardens with decreasing metallicity, due to decreased line blocking 
in the stellar atmosphere winds. This effect is primarily significant above ~ 2.5 Ryd. At 
these energies the blocking is dominated by many Fe IV, Fe V and Fe VI fines, which are 
mainly optically thick near the wind sonic point and optically thin in the outer wind region. 
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As the metallicity is increased, the Fe hnes remain optically thick to higher velocities, thus 
enhancing the blocking influence in this frequency regime. However, below ~ 2.5 Ryd the 
blocking is dominated by C, N, O, Ne, and Ar. The corresponding lines are mainly optically 
thick in the entire wind region. Changes in metallicity therefore do not significantly affect 
the emergent (blocked) fluxes below 2.5 Ryd. Therefore, varying the metallicity has a larger 
effect on the [Ne III]15.6/im/[Ne II]12.8/xm ratio than on the [Ar III]8.99Aim/[Ar II]6.99Aim 
ratio as indicated in Figure 10. For example, for a 40,000 K star increasing the nebular and 
stellar metaUicty from 0.5 to 2 Zq reduces the [Ar III]8.99//m/[Ar Il]6.99/im ratio by a factor 
of ~ 5 and the [Ne IIl]15.6//m/[Ne Il]12.8//m ratio by a factor of ~ 20. 

Figure 10 also demonstrates that the excitation in the LMC nebulae (indicated by star 
symbols) is plausibly enhanced, at least in part, by the reduced metallicity of the LMC. 
However, the low-excitation of the inner Galaxy sources (indicated by the crosses) is prob- 
ably not due to their large metallicities, since these objects appear to have relatively large 
[Ne IIl]15.6/xm/[Ne Il]12.8/xm ratios for their [Ar IIl]8.99/im/[Ar Il]6.99/im ratios, whereas 
increasing the metallicity should selectively reduce the [Ne III]15.6/im/[Ne Il]12.8;um ratio. 
The low excitation of the inner Galaxy sources is most likely due to lower effective tem- 
peratures of the stellar (or cluster) sources. In this respect these objects may resemble the 
typical low-excitation sources observed in the starburst sample of Thornley et al. (2000). 



5. Discussion and Summary 

We have not attempted detailed modeling of individual sources in this paper. However, 
as we have shown, the ensemble of Galactic flne-structure emission line data provide empirical 
support for the generally harder photoionizing SEDs predicted in NLTE wind driven line 
blocked models of hot-star atmospheres. In particular, such models can readily account for 
the presence of high excitation [Ne III] emission hnes in nebular spectra. Earher studies 
which had rehed on LTE atmospheres in the nebular analysis were unable to account for the 
observed [Ne III] emission, and under-predicted the hne intensities by factors of ~ 10 (e.g. 
Baldwin et al. 1991; Rubin et al. 1995; Simpson et al. 1995). Our study largely resolves 
this issue and supports the suggestions of Sellmaier et al. (1996) and Stasinska & Schaerer 
(1997) that this "Ne III problem" was due primarily to a significant under-prediction of Lyc 
photons above ~ 40 eV in the Kurucz LTE atmospheres. 

This conclusion has important consequences for the interpretation of extragalactic fine- 
structure line observations. Thornley et al. (2000) carried out detailed starburst model- 
ing of the [Ne IIl]15.6;um/[Ne Il]12.8;um ratio expected from H II regions ionized by clus- 
ters of stars in which the stellar SEDs were modeled using the new generation NLTE 
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atmospheres. Given that the hottest stars in such models easily produce large nebular 
[Ne IIl]15.6/im/[Ne Il]12.8/im ratios, the low ratios actually observed led to the conclusion 
that the relative number of hot stars is small due to aging of the starburst systems. Our 
analysis supports the conclusion that low [Ne IIl]15.6;um/[Ne Il]12.8/um ratios require the 
removal of the hottest stars as dominant contributors to the ionization of the starburst 
galaxies. 

We summarize our key results. We have compiled archival mid-lR ISO-SWS fine- 
structure lines and H II recombination line observations from H II regions in the Galaxy, the 
LMC, SMC, and M33. We find that: 

1. There is a steep outward increase in excitation in the Galaxy as probed by the fine- 
structure line ratios [Ne III]15.6/im/[Ne II]12.8/im, [Ar III]8.99/xm/[Ar II]6.99/im, and 
[S IIl]18.7//m/[S IIl]33.5//m. Sources in the central kpc of the Galaxy show the lowest 
levels of excitation, and the LMC, and SMC sources show the highest excitation level. 

2. There is a slight decrease in metallicity from ~ 2Zq in the Galactic center to ~ IZq 
in the outer disk. 

3. The increase in nebular excitation across the Galactic disk is likely due mainly to a 
rise in the typical effective temperatures of the exciting stars rather than a hardening 
of the stellar SEDs and enhanced nebular excitation associated with the decline in 
metallicity at fixed effective temperature. 

4. By comparing models to observations, we find that the NLTE wind model atmospheres 
represent an impressive improvement in the realistic description of hot star ionizing 
spectra (compared to LTE static models). In particular the [Ne III] emission observed 
in Galactic H II regions is readily accounted for by these atmospheres. 

5. We estimate an "integrated" [Ne IIl]15.6/xm/[Ne Il]12.8^m ratio of 0.8 for the Galaxy, 
well inside the range of values found for starburst galaxies. The good fit of the NLTE 
models to the data, support the suggestion by Thornley et al. (2000) that the observed 
[Ne III]15.6/im/[Ne II]12.8/im ratio in starburst galaxies is small due to aging effects. 
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List of Figures 

1 Histograms of the observed logarithmic [Ne IIl]15.6//m/[Ne Il]12.8/xm ratio 
(not corrected for extinction) for sources in the our sample - the Galaxy, the 
LMC, the SMC and M33, and also of extra-galactic sources from Thornley 
et al (2000). Filled bars designate sources whose ratios are uncertain. For 
the Milky- Way and the starburst sample, these are all upper limits, for the 
LMC and SMC - lower limits, whereas in M33 both cases are mixed together. 
Our estimation of the integrated [Ne III]15.6/im/[Ne II]12.8/im ratio for the 
Galaxy is over-plotted as vertical long-dashed lines both in the Galactic and 
extra-Galactic histograms. The median [Ne IIl]15.6/xm/[Ne Il]12.8//m ratio 
for the Galaxy is shown as the vertical short-dashed line in the histogram of 
the Galaxy. The Galaxy appears to have higher excitation ratio than most 
starburst galaxies 26 

2 Excitation ratios vs. galactocentric radius for Galactic sources. The plot for 
[Ar IIl]8.99/im/[Ar Il]6.99/im is in the upper panel. In the middle panel we 
show the [Ne III]15.6;um/[Ne II]12.8;umplot, and [S IV]10.5/im/[S III]33.5;um is 
in the lower panel. A rise in the excitation towards the outer Galaxy is clear. 
Solid lines are fits to the gradients (for 2 < Rgai < 11 kpc), whose values 
are given in the text. We designate limits on the given quantities by small 
arrows. The larger arrows indicate lower limits resulting from uncertainties 

in the extinction correction 27 
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3 Top panel - [Ar IIl]8.99//m/[Ar Il]6.99//m vs. [Ne IIl]15.6//m/[Ne Il]12.8//m 

for the entire sample. The data points are represented by the source identifi- 
cations in Tables 2, 5 and 6. A representative error bar for the data is shown 
in the lower right corner. Bottom panel - same as above, with data points 
marked with symbols according to their Galactocentric radius. Pluses desig- 
nate the inner kpc in the Galaxy {Rgai < 1 kpc), circles designate Galactic 



disk sources {Rgai > 1 kpc), and asterisks designate extra-galactic sources. . 28 

4 The Bra/Br/? ratio (upper panel) and the [S IIl]18.7/im/[S IIl]33.5/im ratio 
(lower panel) for Galactic sources in our sample that have measured fluxes for 
these lines, as a function of the galactocentric radius 29 

5 Argon and neon abundances vs. galactocentric radius (left panels) and vs. 
the corresponding excitation ratio (right panels). Solid lines are fits to the 
gradients (for Rgai <18 kpc), whose values are given in the text 30 

6 Comparison between the LTE and NLTE SEDs, for T^ff = 35,000K and 
log (7 = 4.0. The ionization edges of the relevant ions are indicated by vertical 
dashed lines 31 



7 [Ar IIl]8.99/im/[Ar Il]6.99;um vs. [Nc IIl]15.6/xm/[Ne Il]12.8/xm diagnostic 
diagram. Bold faced numbers indicate models. Triangles are 1 Zq NLTE 
models. The dwarfs sequence is designated by solid line, and the super-giants 
sequence is designated by dotted line. Squares are 1 Zq LTE models. Models 
are designated by their effective temperature in lO^K. A representative error 



bar for the data is shown in the lower right corner 32 

8 The fraction of photons that can ionize Ar+, Q argon/ Q vs. the fraction for 
Ne"*", Qneon/Q, for cach one of the model SEDs. The plot follows exactly the 
behavior of the models in the excitation ratios diagram. Figure 7 33 

9 [Ar III]8.99/im/[Ar II]6.99//m vs. [Ne III]15.6yum/[Ne II]12.8;um diagnostic 
diagram with the NLTE dwarfs models (upper panel) and the LTE models 
(lower panel) of various gas densities over-plotted. Triangles are 3,000 cm~^ 
models. Squares are 800 cm~^ models, and diamonds are 100 cm~^ models. 
All model sequences are composed of models with effective temperatures in 

the range 35 (bottom left) to 45 (top right) xlO^K 34 
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10 [Ar IIl]8.99//m/[Ar Il]6.99//m vs. [Ne IIl]15.6//m/[Ne Il]12.8//m diagnostic di- 
agram for dwarfs models with various metallicities. Triangles are solar metal- 
licity models. Diamonds and squares indicate 0.5 Zq and 2 Zq, respectively. 
Models are designated by their effective temperature in lO^K 35 
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Table 1: IR fine-structure lines, which were measured in the spectra selected from the ISO 
archive, with their wavelengths, ionization potentials, and critical densities. 



Line 


Wavelength 


Ionization Potential 


Critical Density 




(/xm) 


(eV) 


(cm~^) 


Bra 


4.05 






Br/3 


2.63 






[Ar II] 


6.99 


27.63 


4.17x10^ 


[Ar III] 


8.99 


40.74 


2.60x10^ 


[Ne II] 


12.81 


40.96 


7.00x10^ 


[Ne III] 


15.56 


63.45 


2.68x10^ 


[Ne III] 


36.01 


63.45 


5.50x10^ 


[S III] 


18.71 


34.79 


2.22x10^ 


[S III] 


33.48 


34.79 


7.04 xlO^ 


[SIV] 


10.51 


47.22 


5.39x10^ 



Table 2: Galactic HII regions data. The fluxes are in units of 10 W cm ^. Asterisks in column 8 mean that there is no unique Rkr, 
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Table 2: Continued. 
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Table 3: Abbreviation for the references of Tables 2. 



1- Afflerbach et al. 1996 

2- Blitz, Fich, & Stark 1982 

3- Brand et al. 1984 

4- Braz & Epchtein 1983 

5- Bronfman, Nyman, & May 1996 

6- Caswell & Haynes 1987 

7- Caswell et al. 1989 

8- Churcliwell, Walmsley, & Cesaroni 1990 

9- Huang et al. 1999 



10- Lang et al. 1999 

11- Lockman 1989 

12- Slieplierd & Churchwell 1996 

13- Slysli et al. 1994 

14- Szymczak, Hrynek, & Kus 2000 

15- Walsh et al. 1997 

16- Wink, Altenhoff, & Mezger 1982 

17- Wink, Wilson, & Bieging 1983 

18- Wouterloot & Brand 1989 



Table 4: Complex regions, whose intensities were summed up. 



OBJECT 


OBSERVATIONS 
(Prom Table 2) 


Arches Cluster 


59-61 


Carina Nebula 


22-25 


IRAS 17431-2846 


62, 65 


IRAS 16172-5134 


41, 42 


IRAS 17424-2859 (Sagittarius A IRS 24) 


55-58 


IRAS 17430-2848 (Pistol star) 


63,64,66 


Omega Nebula (M17) 


73-81 


NGC 6334 


47, 48 


Orion Nebula (M42) 


9-19 


GAL 359.13-00.09 


53. 51 



Table 5: Magellanic Clouds HII regions data. The fluxes are in units of 10 W cm ^. 





NAME 


R.A. (2000) DEC. 


Bra 


Br/3 


[ArII]6.9 


[ArIII]8.9 


[NeII]12.8 


[NeIII]15.6 


[NeIII]36.0 


[SIII]18.7 


[SIII]33.5 


[SIV]10.B 


(1) 


(2) 


(3) 


(4) 


(B) 


(6) 


(7) 




(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


LMC 


113 


30DOR-R136/ 


05 38 42.6 


-69 06 02.5 


2.2 


1.3 








5.2 


27.0 










114 


30DORADUS#2/ 


05 38 35.4 


-69 05 41.2 


5.2 


2.5 


< 0.7 




7.5 


9.0 


64.7 


12.5 


34.7 


98.5 


44.0 


115 


30DORADUS#4.1/ 


05 38 54.2 


-69 05 15.3 


2.8 


1.3 


0.3 




4.9 


6.7 


37.1 


8.7 


22.7 


76.4 


17.9 


116 


30DORADUS#4^/ 


05 38 54.2 


-69 05 15.3 


2.3 


< 1.3 


< 1.8 




4.9 


10.7 


36.0 


12.5 


24.4 


71.9 


1.9 


117 


30DORH2/ 


05 38 44.5 


-69 04 57.8 












15.2 


109.0 










118 


NllA/ 


04 B7 16.2 


-66 23 18.4 


0.3 


0.6 


< 0.4 




1.2 


2.3 


5.1 


0.5 


4.1 


6.3 


2.4 


119 


N1B9-B/ 


05 40 02.4 


-69 44 33.4 


1.2 


0.5 


0.5 




2.0 


6.0 


17.5 


2.2 


8.1 


21.1 


5.1 


120 


N160A-IR/ 


05 39 43.7 


-69 38 30.4 


3.7 


2.5 


< 1.3 




5.8 


16.4 


46.1 


5.9 


21.9 


47.6 


15.8 


121 


N160A1/ 


05 39 43.3 


-69 38 51.4 


4.7 


2.4 


0.5 




7.2 


13.0 


51.0 


7.3 


30.2 


64.8 


25.6 


122 


N160A2/ 


05 39 46.1 


-69 38 36.6 


5.4 


2.8 


1.3 




7.4 


13.9 


38.0 


5.6 


26.5 


49.4 


16.4 


123 


N160A3/ 


05 39 43.7 


-69 38 30.4 


3.7 


1.8 


1.3 




5.1 


12.7 


32.0 


3.7 


18.6 


40.7 


11.5 


124 


N4A/ 


04 52 08.4 


-66 55 23.4 


1.7 


1.0 


1.0 




2.4 


5.0 


13.7 


2.5 


9.4 


21.7 


6.7 


125 


N79A/ 


04 51 47.4 


-69 23 47.7 


< 0.1 


< 0.1 


< 0.4 




< 0.2 


< 0.5 


0.9 


< 0.2 


0.6 


3.3 


0.2 


126 


N83B/ 


04 54 25.2 


-69 10 59.8 


1.0 


0.5 


< 0.4 




1.0 


2.3 


3.5 


0.5 


3.3 


5.5 


1.9 


SMC 


127 


N66/ 


00 59 03.7 


-72 10 39.9 


< 0.2 


0.1 


< 0.6 




<0.4 


< 0.3 


1.3 


< 0.5 


0.7 


1.8 


1.2 


128 


N81/ 


01 09 13.6 


-73 11 41.1 


< 0.6 


0.3 


< 0.3 




0.4 


0.2 


2.7 


0.4 


1.6 


2.4 


1.4 


129 


N88.1/ 


01 24 08.0 


-73 09 02.5 


1.0 


0.8 


< 0.6 




0.8 


< 0.4 


3.0 


0.5 


1.1 


0.9 


4.9 


130 


N88.2/ 


01 24 08.1 


-73 09 03.1 


1.5 


< 2.9 


< 1.5 




<1.2 


< 1.6 


6.7 


< 3.5 


< 3.3 


< 3.5 


6.8 



Table 6: M33 Clouds HII regions data. Rm33 is the M33 angular galactocentric radius in arc minutes. The fluxes are in units of 10 W cm 





NAME 




R.A. (2000) DEC. 


^M33 


Bra 


Br/3 


[Aril] 6. 9 


[Arin]8.9 


[NeII]12.8 


[NeIII]15.6 


[NeIII]36.0 


[SIII]18.7 


[SIII]33.5 


[Siv]ia.5 


(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


131 


M33C004/ 


01 


33 59.3 


+30 35 49.2 


3'21" 




< 0.2 






0.5 


< 0.6 










132 


M33C027/ 


01 


33 45. .5 


+30 36 51.1 


2'22" 




< 0.2 






0.2 


< 0.2 










133 


M33C041/ 


01 


33 35.5 


+30 39 30.5 


2'22" 




< 0.1 






0.2 


< 0.5 










134 


M33C045/ 


01 


33 28.8 


+30 40 24.7 


3'45" 




< 0.2 






0.5 


1.1 










135 


M33C049/ 


01 


33 33.9 


+30 41 28.3 


3'45" 




< 0.2 






1.2 


0.8 










136 


M33C079/ 


01 


34 00.2 


+30 40 48.6 


1'41" 




< 0.2 






0.8 


0.6 










137 


M33C088W/ 


01 


34 15.5 


+30 37 12.4 


3'21" 




< 0.2 






0.2 


0.3 










138 


M33C095/ 


01 


34 10.9 


+30 36 18.1 


3'45" 




< 0.2 






0.3 


< 0.4 










139 


M33C214/ 


01 


33 29.1 


+30 31 46.8 


7'52" 




< 0.2 






0.4 


< 0.2 










140 


M33C2B1/ 


01 


33 36.6 


+30 20 13.4 


18\55" 




< 0.2 






< 0.2 


0.3 










141 


M33C280/ 


01 


32 4B.2 


+30 38 53.7 


8'48" 




< 0.2 






< 0.2 


0.6 










142 


M33C302/ 


01 


34 06.9 


+30 47 25.9 


8'38" 




< 0.2 






0.6 


< 0.5 










143 


M33C623.1/ 


01 


33 16.5 


+30 52 49.1 


14'38" 




< 0.2 






0.3 


0.9 










144 


M33C623.2/ 


01 


33 16.3 


+30 52 47.8 


14'38" 


0.1 


< 0.1 


< 0.2 


0.1 


< 0.4 


0.7 


< 0.7 


0.3 


< 0.3 


0.4 


145 


M33C638E/ 


01 


33 16.3 


+30 56 43.8 


18'23" 




< 0.2 






< 0.2 


0.3 










146 


M33C680B/ 


01 


34 33.5 


+30 46 50.5 


9'21" 




< 0.3 






1.8 


1.7 










147 


M33C680C/ 


01 


34 32.1 


+30 46 59.9 


9'21" 




< 0.2 






2.5 


2.7 










148 


M33C691/ 


01 


34 16.4 


+30 51 55.1 


13'13" 




< 0.3 






0.2 


0.8 










149 


M33C710/ 


01 


34 13.6 


+30 33 42.5 


5'49" 




< 0.2 






0.4 


0.2 











-42 - 



Table 7: The model atmospheres parameters - stellar effective temperature, surface gravity, 
stellar radius, wind parameters (M and Voo) and hydrogen photoionizion rates of NLTE 
typical dwarfs. The parameters are given for half-solar, solar, and twice-solar metallicities. 











0.5Zq 












2Zq 






log 9 


R* 


logQ 


Vao 


lO-^'M 


logQ 


foe 


lO-'^M 


logQ 




lO-'^M 


(K) 


(cm s~^) 


(Rq) 


(s-l) 




(M0j/r-l) 






{M02/r-i) 






(M0yr-l) 


35,000 


4.0 


9.0 


48.37 


2000 


0.8 


48.35 


2004 


1.6 


48.37 


2008 


3.4 


36,000 


4.0 


9.2 


48.52 


2079 


1.0 


48.51 


2093 


1.9 


48.53 


2086 


3.9 


37,000 


4.0 


9.4 


48.67 


2158 


1.2 


48.68 


2177 


2.3 


48.68 


2181 


4.6 


38,000 


4.0 


9.6 


48.76 


2283 


1.5 


48.73 


2272 


2.8 


48.73 


2277 


5.2 


39,000 


4.0 


9.8 


48.88 


2358 


1.8 


48.87 


2362 


3.3 


48.85 


2365 


6.1 


40,000 


4.0 


10.0 


48.99 


2450 


2.2 


48.98 


2450 


4.0 


48.98 


2460 


7.0 


41,000 


4.0 


10.2 


49.09 


2548 


2.8 


49.08 


2541 


4.8 


49.07 


2564 


8.3 


42,000 


4.0 


10.4 


49.19 


2637 


3.4 


49.18 


2640 


5.7 


49.16 


2620 


9.6 


43,000 


4.0 


10.6 


49.26 


2732 


4.0 


49.26 


2719 


6.7 


49.26 


2722 


11.0 


44,000 


4.0 


10.8 


49.34 


2803 


4.9 


49.33 


2806 


7.8 


49.32 


2800 


13.0 


45,000 


4.0 


11.0 


49.42 


2894 


5.8 


49.40 


2923 


9.1 


49.39 


2900 


14.0 



Table 8: The model atmospheres parameters of our NLTE typical giants for half-solar, solar, 
and twice-solar metallicities. 











O.5Z0 






IZq 






2Zo 






logs 


R* 


logQ 


foo 


lO-^'M 


logQ 


Voo 


lO-^'M 


logQ 




lO-'^M 


(K) 


(cm 


(Re) 


(s-l) 


(fcms~^) 


(M0:yr-l) 


(s-i) 




(M0?;r-l) 


(s-l) 




(M0yr-l) 


35,000 


3.25 


24.0 


49.54 


1988 


48 


49.54 


1998 


63 


49.52 


2008 


83 


36,000 


3.30 


23.4 


49.60 


2058 


52 


49.60 


2061 


69 


49.59 


2083 


91 


37,000 


3.35 


22.8 


49.64 


2122 


54 


49.63 


2125 


74 


49.62 


2128 


99 


38,000 


3.40 


22.2 


49.65 


2172 


58 


49.67 


2180 


78 


49.66 


2193 


110 


39,000 


3.45 


21.6 


49.72 


2227 


61 


49.69 


2235 


83 


49.70 


2240 


110 


40,000 


3.50 


21.0 


49.74 


2313 


64 


49.75 


2311 


88 


49.74 


2329 


120 


41,000 


3.55 


20.4 


49.79 


2364 


66 


49.76 


2367 


91 


49.77 


2376 


130 


42,000 


3.60 


19.8 


49.81 


2424 


68 


49.78 


2426 


95 


49.79 


2431 


140 


43,000 


3.65 


19.2 


49.84 


2493 


70 


49.85 


2483 


100 


49.82 


2490 


140 


44,000 


3.70 


18.6 


49.86 


2556 


71 


49.86 


2548 


100 


49.85 


2569 


150 


45,000 


3.75 


18.0 


49.87 


2609 


73 


49.88 


2610 


110 


49.87 


2618 


160 



-43- 



Table 9: The model atmospheres parameters of the LTE dwarfs for solar metaUicity. 





log 3 


R* 


logQ 


(K) 


(cm s^'^) 


(Rq) 


(s-i) 


35,000 


4.0 


9.0 


48.43 


36,000 


4.0 


9.2 


48.57 


37,000 


4.0 


9.4 


48.70 


38,000 


4.0 


9.6 


48.81 


39,000 


4.0 


9.8 


48.92 


40,000 


4.5 


10.0 


48.96 


41,000 


4.5 


10.2 


49.05 


42,000 


4.5 


10.4 


49.14 


43,000 


4.5 


10.6 


49.23 


44,000 


4.5 


10.8 


49.31 


45,000 


4.5 


11.0 


49.38 



